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Abstract 

Multi-stage sensing is a novel concept that refers to a general class of spectrum sensing algorithms 
that divide the sensing process into a number of sequential stages. The number of sensing stages and 
the sensing technique per stage can be used to optimize performance with respect to secondary user 
throughput and the collision probability between primary and secondary users. So far, the impact of 
multi-stage sensing on network throughput and collision probability for a realistic network model is 
relatively unexplored. Therefore, we present the first analytical framework which enables performance 
evaluation of different multi-channel multi-stage spectrum sensing algorithms for Opportunistic Spectrum 
Access networks. The contribution of our work lies in studying the effect of the following parameters on 
performance: number of sensing stages, physical layer sensing techniques and durations per each stage, 
single and parallel channel sensing and access, number of available channels, primary and secondary 
user traffic, buffering of incoming secondary user traffic, as well as MAC layer sensing algorithms. 
Analyzed performance metrics include the average secondary user throughput and the average collision 
probability between primary and secondary users. Our results show that when the probability of primary 
user mis-detection is constrained, the performance of multi-stage sensing is, in most cases, superior 
to the single stage sensing counterpart. Besides, prolonged channel observation at the first stage of 
sensing decreases the collision probability considerably, while keeping the throughput at an acceptable 
level. Finally, in realistic primary user traffic scenarios, using two stages of sensing provides a good 
balance between secondary users throughput and collision probability while meeting successful detection 
constraints subjected by Opportunistic Spectrum Access communication. 
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I. Introduction 

Opportunistic Spectrum Access (OSA) radios are wireless communication devices that exploit tem- 
porarily vacant licensed frequency bands 121, [Sl- To be able to use these spectral vacancies, OSA radios 
need to sense the activity of the licensed user in the frequency bands of interest. As spectrum sensing is 
the enabling technology for OSA, many sensing techniques have been proposed in the literature covering 
a variety of design aspects H. This includes different sensing mechanisms at the physical layer as well 
as various collaborative sensing techniques and measurement fusion methods at the link layer. 

In a conventional spectrum sensing algorithm, a OSA radio (denoted later as Secondary User (SU) 
radio) has to periodically scan the licensed spectrum for any incoming Primary (licensed) Users (PUs), 
and detect the PUs with a predetermined probability. When the SU radio's spectrum sensing physical 
layer triggers an alarm that a PU exists, i.e. due to a successful detection or a false alarm, the SU radio 
should stop transmission on the respective channel. However, the sensing procedure might result in a 
false alarm, thus, reducing the channel access opportunity for the SU radio. The probability of a false 
alarm can be decreased by, e.g. increasing the sensing time, however, this might lead to a decrease in 
the throughput achieved by the SU device, as demonstrated in e.g. m-Q. Besides, the probability of 
false alarm can be decreased by altering the decision threshold regarding the presence of the primary 
user. However, this leads to an increase in the probability of a collision between the SU radios and PU 
devices. This in turn may lead to a decease in throughput caused by packet losses at the SU radio. 

One of the ways to increase channel utilization in OSA communication is to use multi-stage sensing. In 
multi-stage sensing, an alarm (due to successful detection or false alarm) on the PU presence is followed 
by 5 — 1 cycles of sensing (and potential transmission). Thus, the presence of a PU would be declared, 
and the SU radio would vacate a channel, if and only if S consecutive alarms are generated. The sensing 
technique used at each of the S stages can be independent from one another. Accordingly, multi-stage 
sensing provides additional degrees of freedom to optimize performance metrics, namely, throughput of 
SU radio and the collision probability between PUs and SUs. 

A. Related Work 

Multi-stage sensing has been introduced formally in the IEEE 802.22 standard lUl, ||9l, with 5 = 2, 
while the first analyses of a multi-stage sensing algorithm have been performed in |[TOl . ifTTI (with 
S* = 2 in case of [111 ) only for a single PU channel. Note that ifTOll . ifTTI did not consider the random 
nature of the OSA traffic which has a notable effect on the performance of the OSA network. Authors 
in |[T6ll extended their previous work ifTOl to the case of two physical layer sensing techniques: (i) energy 
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detection and (ii) feature (cyclostationary) detection. Furthermore, perfect detection at the last stage of 
sensing and zero delay in channel switching is not assumed. In their model, the frequency of sensing 
the channel increases when the sensing procedure indicates that PU is possibly present. Finally, in | fT2l | 
a multi-channel sensing technique with S = 2 and one sensing algorithm was considered, where each 
sensing phase occurred consecutively in the same time sloj^. 

B. Our Contribution 

We conclude that the literature lacks a unified analytical framework that evaluates multi-stage spectrum 
sensing options when S > 2, considering wide range of algorithms and parameter setups, sensing 
techniques, various traffic conditions for SU and PU radios, SU buffer availability and multi-channel 
access. In this paper we provide a comprehensive analytical framework that models the multi-stage 
multi-channel spectrum sensing concept. 

Our work is specifically related to the recent work of |fT6l, but should be treated as a new approach for 
analyzing multi-stage spectrum sensing algorithms. In contrary to |16|, we provide a more realistic model 
for PU activity, where we consider the PU traffic, not only in the channel where the SU is operating, 
but also in all available channels. Furthermore, we account for the scarcity in the bandwidth available 
for OSA operation where our model accounts for a limited pool of available PU channels. Moreover, 
our model applies to OSA devices that have a single radio operating on a narrowband channel and OSA 
devices that are defined to have a stack of narrowband radios that can utilize or sense multiple PU 
channels at the same time. In contrast to |[T6l , our model accounts for different traffic models for the 
OSA network, including constant and variable bit rate streams. Besides, the model covers unbuffered and 
buffered OSA systems. Performance is quantified by throughput and a new metric, that is not considered 
in |[T6l . which is the collision probability between the PU traffic and SU devices. We use the proposed 
analytical framework to evaluate three new multi-stage sensing algorithms that aim at optimizing the 
aforementioned metrics. In the context of these algorithms, we introduce a new sensing mode, denoted 

'Note that the idea of multi-stage sensing is conceptually similar to the multi-dwell algorithms used in CDMA code acquisition 
analyzed in 1131 Sec. II] or in UWB phase tracking analyzed in 1141 Sec. II], II5I Sec. II]. However, there are two main differences 
between multi-stage sensing and multi-dwell acquisition that require a novel modeling approach for multi-stage spectrum sensing. 
Firstly, in the multi-dwell algorithm, sensing is performed at each stage in parallel, i.e. multiple components sense the same signal 
at the same time. This is contrary to multi-stage spectrum sensing where parallel detection is problematic due to the hardware 
cost. Secondly, the detected signal for multi-dwell acquisition is stationary in the time domain, which cannot be guaranteed in 
OSA systems, due to the traffic fluctuation of the PU. 
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as pre-sensing, which targets minimizing the coUision probability. Given that |[T6l already evaluated the 
effect of using different physical sensing methods, and channel switching delay, we do not include it in 
our study. Since our introduced model accounts for the randomness of the PU and SU traffic, and the 
most widely used method of analyzing network traffic is via Markov chain analysis [17, Ch. 11], we 
follow the same approach here, just like IfTOl . |[T6l . 

The rest of the paper is organized as follows. The system model is introduced in Section [III The 
analytical model is presented in Section |lIIJ with the multi-stage sensing algorithm design examples 
analyzed in Section |IVl The numerical results are presented in Section |Vl Finally, the paper is concluded 
in Section IVll 

II. System Model 

A. General Assumptions 

SUs communicate opportunistically over multiple narrowband channels of throughput W^kbps, that 
are randomly occupied by PUs (all variables are summarized in Table We consider an OSA network 
composed of two users where one is always the transmitter and the other one is always the receiver. The 
reduction of the OSA network to only two nodes allows us to neglect the impact of collisions of random 
access control packets on the performance of spectrum sensing algorithms. Note that the same assumption 
has been used in an earlier study of multi-stage sensing ifTOl . Our model can be extended to the multiple 
nodes case, however, for the brevity of the discussion, this is not addressed here. Furthermore, PUs and 
SUs, are assumed to arrive and depart at discrete times that are multiples of T, where T is the slot length 
in seconds. We assume that the transmitting SU transmits the sensing results to the receiving SU through 
the dedicated control channel. Hence, the SU pair is perfectly synchronized. These two assumptions, 
vital in simplifying the analysis, are common in spectrum sensing literature, refer to ifTSl for a good 
overview of related work in this area. Finally, we assume that the time taken by the SUs to switch from 
communicating mode to sensing mode, or from one channel to another, is negligible. 

In this work we focus on OSA link throughput and PU/SU collision analysis. We assume that the 
obtained OSA network throughput is proportional to the duration of time the SUs use the channel and 
the number of channels utilized. Accordingly, throughput decreases if the SU is idle or senses the channel. 
Besides, in cases where the SU and the PU use the same channel simultaneously, we assume that the 
SU frame cannot be decoded and is considered lost. Furthermore, no acknowledgments are considered 
for the data transmitted by SUs. Finally, the frame error rate is assumed to be zero on PU free channels. 
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TABLE I 

Summary of Variables Used in the Paper 



Variable 


Description 


Unit 


S, M, N 


number of: sensing stages, sensing radios, available channels 


— 


B, b 


buffer size, buffer utilization 


— 


R, G 


average: SU throughput, collision probability 


bps, — 


W 


channel throughput 


kbps 


3, J 


SU state vector: single narrowband radio case, parallel narrowband radio case 


— 


r, Ts {Va) 


set of SU states: feasible, where the SU is sensing (active) 


— 


e, * 


set of states: single narrowband radio, parallel narrowband radio 


— 


Oi (92) 


set of states for single narrowband radio: SU sensing, PU exists (PU absent) on channel 


— 


6'v, '(Aw 


state: single narrowband radio, parallel narrowband radio 


— 


T, Ts, Tt 


slot length, sensing duration at stage S, sensing duration for quiet and pre-sensing 


(m)s, (yu)s, (/j)s 


p, s 


variables denoting: PU, SU 


— 


i, X, k, I 


supporting variables 


— 


c, m 


indices: channel, radio 


— 


ria, n-a 


number of chaimels with PU: arrivals, no arrivals 


— 




number of channels with PU: departures, no departures 


— 


Pf,s (Pm,s) 


probability of false alarm (mis-detection) at sensing stage 


— 


Pf,t (Pm,t) 


probability of false alarm (mis-detection) at pre-sensing and quiet modes 


— 


Px,a (Px.d) 


probability of arrival (departure) for user x 


— 


Wi, Vi 


vector of variables at slot i: single narrowband radio, parallel narrowband radio 






supporting vector of variables 




X, y 


variables for conditional probabiUties: single narrowband radio, parallel narrowband radio 




f 


SU traffic state for single channel SUs 




F 


vector of SU traffic state for parallel narrowband radio 




I, T 


vector of PU states, set of all permutations of / 




TT, A 


vector of stationary probabilities, transition probability matrix 




Fn^ 


number of frames generated at slot i — 1 for parallel narrowband radio 




Fxi, Fsi 


total number of frames at slot i: available for transmission, buffered 






total number of radios at slot i: available for transmission, active 




f^Aj, ^iQi 


set of radios at slot i: active, idle, quiet 
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B. SU and PU Traffic Model 

We assume that the SU sends packet traffic randomly (either with a constant or a variable bit rate), 
where packets are divided into frames. In the idle mode, the SU has no frames to transfer, while in the 
active mode, the SU is either transmitting frames or sensing the channels. The probabilities of arrival 
and departure of a SU frame are denoted by ps^a and ps^d, respectivelj^. Further, we consider SU nodes 
with and without data buffering capabilities. In the case of unbuffered SU nodes, incoming frames are 
discarded when no channels are available for transmission. In a buffered node case, a buffer of size B 
frames is present that stores frames only when the node is in a mode of operation where sensing occupies 
the whole time slot. New frames are discarded if the buffer is full. The presence of the buffer increases 
the average throughput as incoming frames are not discarded when SUs are sensing the channel for the 
whole slot length. 

PU channel occupancy is also time varying. That is, the probabilities of the start and end of a channel 
occupancy are denoted by pp^a and pp rf, respectively. Finally, the average PU occupancy duration and 
idle times are assumed to be the same for every PU channel. 

C. Multistage Sensing Algorithms 

When a SU transmitter generates a new packet, it first sends a request to send the packet on the 
dedicated control channel to the SU receiver. Both the sender and the receiver are constantly tuned to 
the control channel during the idle mode. The connection establishment time is assumed to be small in 
comparison to the subsequent transmission time and can therefore be neglected. 

Once connection is established, the sending SU node starts with the first stage of the multi-stage 
sensing algorithm on the first PU channel from a list of available PU channels. In every stage, the SU 
transmitter first senses the channel for Ts < T seconds (except in the special cases, explained later in this 
section, where sensing occupies the whole time slot). Immediately after sensing, in the same slot, a SU 
frame of duration T — Ts seconds is sent regardless of the sensing procedure outcome^. If the outcome 

^Note that the arrival and departure rates of the SU frames are constant and are independent of the PU state and the SU sensing 
mechanism performance throughout the model. This assumption is typical for system level analysis of OSA networks. Inclusion 
of SU traffic which is dependent on the PU state in the model will make the observation of the SU network performance difficult, 
e.g. it will be hard to infer whether the decrease in throughput is a result of decreasing SU traffic or the effect of multi-stage 
sensing algorithm. While the inclusion of the PU-dependent SU traffic in our model is an interesting topic to consider, we leave 
it for future studies. 

^Note that contrary to 1101 . 1161 we do not assume a varying sensing time at each stage. Varying sensing time would complicate 
the packet forwarding process as frames always have a fixed length and cannot adapt to varying sensing duration. 
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of the sensing procedure at stage j indicates that a PU is possibly present, the SU nodes proceed to the 
next sensing stage. If the outcome of the sensing procedure at stage j indicates that the PU is possibly 
absent, the SUs start a new multi-stage sensing cycle. 

The probabilities of false alarm and mis-detection for any of the S stages of multi-stage sensing are 
denoted by pf ,, and Pm,s, respectively. Note that, generally, increasing Ts decreases pf^g and Pm,s- On 
the other hand, this leads to a decrease in the SUs throughput, as more time is spent on sensing for PU 
activity. Also note that SUs can transmit during each stage of the sensing cycle even though the PU is 
present as a result of a mis-detection, which is also the case for single-stage sensing algorithms. 

The multi-stage spectrum sensing algorithm might outperform the single-stage sensing counterpart as a 
result of increased false alarm probability. In single-stage sensing algorithms, a single false alarm forces 
the SU to stop transmission immediately, while in the case of multi-stage spectrum sensing, the SU 
checks in S consecutive sensing stages if a PU is actually on the channel without stopping the on-going 
communication. This might improve the SUs' utilization of channel vacancies specially in scenarios where 
the false alarm probability is high and/or the PU traffic is slow. This is further explored in Section jV] 

In the multi-stage spectrum sensing algorithm, the SU transmission continues on the respective channel 
until the SU detects the PU presence in all S stages of sensing till the last stage. The SUs may then 
switch to the next channel and restart the multi-stage spectrum sensing algorithm, or stay on the same 
channel in a prolonged sensing stage that involves no transmission, depending on the radio architecture 
as explained in the following sections. 

I ) Multi-channel Sensing using a Single Narrowband Radio: First, we focus on OS A nodes with only 
one radio for sensing or communication which can operate only on one frequency band at a time. SUs 
can switch to a different narrow frequency band when they decide to vacate their current band due to 
the possible presence of a PU. 

The proposed model serves as a general framework enabling the analysis of any multi-stage sensing 
algorithm. Due to the infinite number of possibilities of multi-stage sensing algorithms, we focus our 
analysis on the algorithm presented in the IEEE 802.22 standard 13, as well as our proposed extensions. 
The basic algorithm operates as follows. When the SUs reach the last stage of multi-stage sensing and the 
outcome of the sensing procedure indicates that a PU is possibly present, the SUs switch to a prolonged 
channel observation stage. We denote this stage by the quiet mode. In the quiet mode, the SUs stay on that 
channel for the whole slot duration Tt = T, to sense for PUs. If the outcome of the sensing procedure at 
the quiet mode indicates that a PU is possibly present, the SUs switch to another PU channel and begin 
a new multi-stage sensing cycle. Otherwise, the SUs stay on the same channel and restart the multi-stage 
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sensing cycle to detect incoming PUs or a PU that was mis-detected. Regarding the quiet mode, a false 
alarm with probability pj^t, or a mis-detection with probability pm,t can happen (just as in [16], but 
contrary to ffTOl , where detection at this mode of operation is assumed perfect). Since Tf > Tg, then pj^ 
and p„i,t can be chosen to be less than pj^g and Prn,s, respectively. This means that the sensing procedure 
for the prolonged channel observation stage, also called fine sensing in the IEEE 802.22 standard [8], has 
higher precision than that for stages 1 to S. However, this comes at the expense of a potential decrease 
in channel utilization by the SU. The motivation behind the quiet mode is to decrease the probability of 
the SU leaving a vacant channel based on false alarms (the impact of the quiet mode on the multi-stage 
spectrum sensing performance is discussed in more detail in Section jV]). 

In addition to the aforementioned algorithm, we propose and analyze algorithms which involve the 
usage of a prolonged channel observation period at the beginning of the sensing cycle. We denote this 
prolonged sensing period by the pre-sensing mode. The length of this period is Tt = T, exactly as in the 
quiet mode. The motivation behind the pre-sensing mode is to decrease the collision probability with the 
PU, since p^^t is small, when the SUs switch to a new channel. At the end of pre-sensing mode, the SUs 
enter the first stage of sensing (if no PU has been detected or a PU has been mis-detected) or switches 
to a new channel and starts a new sensing procedure with pre-sensing (if a PU has been detected or 
a false alarm occurred). The quiet and pre-sensing stages have a considerable effect on throughput and 
collisions as our results demonstrate. Accordingly, we introduce and analyze algorithms combining all 
permutations of these modes of operation. 

2 ) Multi-channel Sensing using Parallel Narrowband Radios: Apart from the single narrowband radio 
case, we consider a system, where both SUs have M narrowband radios that can operate in parallel 
independently. That is, at any given time, each narrowband radio can be idle, sensing or communicating 
regardless of the state of the other radios. Besides, the radios may operate on non-contiguous frequency 
bands. In this model we assume that the number of available radios is equal to the number of available 
channels, N, which is the only relevant case in practice. The incoming frames are divided across the 
radios with no prioritization. Each radio follows the basic single channel multi-stage sensing algorithm 
described in the earlier section, except once the quiet mode is reached, the radios stay in the quiet stage 
until the outcome of the sensing procedure indicates that the channel is vacant. If, in the quiet mode, the 
channel is correctly detected to be vacant or a PU occupying the channel is mis-detected, the sensing 
procedure outcome would indicate that the channel is vacant and the corresponding radio proceeds to the 
first stage of the spectrum sensing algorithm. 
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III. Analytical Model 

We present the foundation for the analytical model for the single and parallel narrowband radio cases 
separately. We divide the discussion into the Markov chain state definitions, state transition probabilities, 
stationary probability calculation and metrics calculation. 

A. Multi-channel Sensing using a Single Narrowband Radio 

1) State and State Transition Definitions: Let 9^ denote the state of the Markov chain, where v = 
{I, f,j,b,c}. In vector v, c G {1, • • • , N} is the operating channel index, which denotes the channel 
the SU is operating on, or was operating on before proceeding to the idle mode. SU traffic status is 
represented by /, where / = 1 when the SU has a new frame to send, and / = otherwise. The SU 
mode of operation is represented by j € F, where F is defined differently for each single narrowband radio 
sensing algorithm. For the proposed algorithms, j = implies that the SU is idle and j € {!,••• , S} 
indicates that the SU is in sensing stage j of multi-stage sensing, while j = S + 1 and j = 5 + 2 are 
the examples of the special cases that in the context of this work indicate that the SU is in the quiet 
mode and the pre-sensing mode, respectively. Define F^ as a subset of F representing the modes of 
operation where the SU is performing multi-stage sensing. Define F^ as a subset of F representing the 
modes of operation where the SU is active, that is, the SU is either performing multi-stage sensing or 
is in the pre-sensing or quiet mode. F^ = {1, • • • , 5"} for all algorithms, however, Ta depends on the 
specific algorithm. The buffer level at the SU is indicated by & € {0, • • • , B}. Finally, the PU status are 
described by /, |/| = N, where I{x) = 1 when a PU operates on channel x, and I{x) = otherwise. 
When transitioning from time slot i to time slot i + 1, the state of the Markov chain at slot i + 1 (the 
SU mode of operation, buffer level and operating channel in specific) depends on the outcome of the 
sensing procedure that tests the existence of a PU during time slot i (not time slot i + 1). Hence, for 
mathematical convenience, / is defined to describe the PU status in the previous time slot and not the 
current time slot. 

The transition probability from state to state 6^^, where Vx = {Ix,fx,jx,bx,Cx} and subscript 
X = {1, 2} for c, /, /, j, b, denotes current and next time slots, respectively, is expressed as 

3 

Fv{X) = U{X)llFvk{X), (1) 

k=l 

where again, for notational convenience, X = {6'v, I 6lvJ. In O, Pri(-) denotes the probability of PU 
status change, Pr2(-) denotes the probability of SU traffic status change, Pr3(-) denotes the probabilities 
of the sensing procedures outcomes, and U{x) is the transition feasibility function and equals to one if 
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X belongs to a specified set of feasible state transitions and U{x) = otherwise. Since the traffic of 
the PU and SU is independent of the sensing algorithm, Pri (•) and Pr2 (•) are the same for all parallel 
narrowband radio algorithms, described later in Section UlI-BI However, as Pr3 (•) and U {■) describe 
transitions in the radios mode of operation, they are unique to the sensing algorithm and are described in 
Section ITV-AI The transition probabilities are consequently used in calculating the stationary probabilities 
as described in the following section. 

2) Stationary Probability Calculation: Define as the set of all feasible Markov chain states. Then, 
define A as the transition probability matrix where A {k, I) denotes the transition probability from state 
O (k) to state @ (/), where /c, / G {1, • • • , |©|}. © is evaluated by creating a vector of unique states, 9^,, 
for all feasible combinations of c,j, b, I, f and directly mapping them to the values of ©. Further, define 
vr as the vector of stationary probabilities where tt {k) is the stationary probability of state Q{k). Then 
vr can be evaluated by solving vr = vrA, knowing that X]l=i ^ (^) = 1- The stationary probabilities are 
then used to calculate the desired metrics. 

3) Performance Metrics Calculation: Based on the stationary probability vector, vr, we compute the 
average throughput obtained by the secondary network, R, and the expected number of collisions between 
the transmitting SU and PUs per time slot, G. Define ©i as the set of states ^v, with j G Ts and / (c) = 1, 
i.e. ©1 is the set of states where the SU is performing multi-stage sensing and a PU was present on the 
channel on which the SU was operating during the previous time slot. Define ©2 as the set of states 6^, 
with j G Ts and / (c) = 0. In calculating R, we assume that the SU frame cannot be decoded in case 
of a collision between the SU and a PU, and accordingly, R can be expressed as 

R = w'^^{ T^{i)Pp4+ J2 vr(i)(l-pp,a) I , (2) 

\0(j)Gei e(j)ee2 / 

where (T — Ts)/T is the portion of a time slot SU spent in transmission. Besides, G can be expressed 
as 

G= ^ 7r(i) (1 + ^ 7r{i)pp^a- (3) 
6(1)661 6(i)e62 

Note that in Q and (|3]l, the information regarding the PU status according to 9^, given by /, pertains to 
that of the previous time slot. Accordingly, the probability that no PUs exist at the current time slot at 
channel c given that /(c) = 1 equals pp^d, i-c the probability that the PU user that was utilizing channel 
c in the previous time slot is not active during the current time slot. Likewise, the probability that a PU 
is active on channel c given that /(c) = 1 is (1 — Pp^d)- This accounts for the Pp^d and (1 — pp^a) terms 
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in ([2]) and the ^ and (1 — Pp^d) terms in Q. In the following sections we will describe probabilities 
governing ([T]) for each sensing algorithm. 

4) Definitions ofPri{-) and Pr2(-)." Let denote the number of channels with PU arrivals, denote 
the number of channels with no PU arrivals, Ud denote the number of channels with PU departures and 
denote the number of channels that remained occupied by the PU. It follows that = /i o I2, 
ria = Ii o I2, rid = h o I2, and n-^ = Ii o I2, where x is the bitwise binary complement of x and o is 
the base- 10 dot product operator. Accordingly, Pri(-) is expressed as 

Pri {X) = p;:ail - Pp,ar%:dil - V,4T\ (4) 
where Ii,/2 S T and T is the set of all 2^ permutations of /. Pr2(-) is expressed as 



Pr2 (A") 



1 -Ps,a, 


fl 


= /2 = 


0, 


Ps,ai 


fl 


= 0,/2 


= 1 




fl 


= l,/2 


= 


1 -Ps,d, 


fl 


= /2 = 


1. 



(5) 



B. Multi-channel Sensing using Parallel Narrowband Radios 

1 ) State and State Transition Definitions: Let ip^ denote the states of the Markov chain, where w = 
{/, F, J, b}. In w^, / and b are defined the same as for the single narrowband radio case. Let tp^. denote 
the Markov chain state at time slot i, where Wi = {Ii, Fi, Ji,bi}. In the context of our work, the 
parallel channel SU is assumed to have M single channel radios, the radio states are described by J, 
|J| = M, where J (x) S P describes the mode of operation of radio x. For the considered algorithm 
r = {0, • • • , S" + 1} and J (x) = S + 1 indicates that radio x is in the quiet mode. The rest of the radio 
modes of operation are defined as for the single narrowband radio case. Furthermore, Pg = {1, • • • , S}, 
as in the case of singe parallel radios, defines the subset of P that represents the modes of operation 
where the SU radio is performing multi-stage sensing. 

SU traffic is assumed to be in the form of a frame stream with frame arrival and departure probabilities 
as in the single narrowband radio case. The difference is that up to M frames can be generated per time 
slot. Thus, for analysis, a slot can be hypothetically segmented to M parts where a new frame can be 
generated at each part. The state of these frames is represented by Fi, \Fi\ = M, where Fi (x) = 1 if 
a frame was created at division x of slot i — 1 and Fi {x) = otherwise. Hence, the total number of 
frames generated at slot i — I, denoted by F^^, can be expressed as F^r^ = X]x=i (^)- These frames 
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are accumulated with buffered frames from time slot i — 1 and are ready for transmission at time slot 
i. The excess frames that cannot be transmitted at time slot i due to the limited number of available 
channels are buffered where a maximum of B frames can be stored. The total number of frames available 
for transmission at time slot i, denoted by Ft^, can be expressed as Ft^ = + Fjv^ where is the 
number of buffered frames from time slot i — I. Radios that are in the quiet mode at slot i are not used 
for frame transmission. Hence, the total number of radios available for frame transmission at slot i can be 
expressed as Mp^ = X]m=i ^ i^i ("^)) where V {Ji (m)) = 1 when Jj (m) / 5 + 1 and V {Ji (m)) = 0, 
otherwise. Accordingly, the total number of radios that will be active at slot i, which equals the total 
number of frames that will be transmitted, is Ma^ = minjFj^^, Mi?.}. These radios are selected from the 
set of radios that are available for frame transmission at slot i in ascending order of the radio indices 
(due to the symmetry of the system). 

Denote the set of radios that will be active at slot i by Q.Ai where 1^2 a J = Ma,- If Mp, > Ma^, the 
extra radios that are available for frame transmission at slot i will be idle. Denote the set of these radios 
by 0.1- where [0/^ [ = max{0, Mp. — Fp^ }. Radios that are in the quiet mode will stay in the quiet mode 
until the outcome of the sensing procedure indicates the absence of a PU. Denote the set of radios that 
are in the quiet mode at slot i by Q.q^ where I^Jq. | = M — Mp.. If Mp^ > Fp,, then Ma, = Fp. and 
all available frames will be transmitted, thus, 6, = 0. Otherwise, M^. < Fp^ and only Ma, frames will 
be transmitted, hence, hi = Tam.{B,Fp. — Ma,}- The rest of the frames, if any, will be dropped. The 
transition probabiUty from state ipv/i to state ip^^ is given by 

3 

Y>v{y) = u{y)\{vvk{y), (6) 

k=l 

where, for notational convenience, y = {V'wal^wi}- In ©, Pri (•) represents the probability of PU status 
change, Pr2 (•) denotes the probability of SU traffic status change, Pra (•) represents the probabilities 
of the sensing procedures outcomes, and U {■) is the transition feasibility function and U{x) = 1 if x 
belongs to a specified set of feasible state transitions and U{x) = otherwise. Since the traffic of the 
PU and SU is independent of the sensing algorithm, Pri (•) and Pr2 (•) are the same for all algorithms. 
However, as Prs (•) and U {■) describe transitions in the radios mode of operation, they are unique to the 
sensing algorithm and are described in Section IIV-BI 

2) Stationary Probability Calculation: Denote as the set of all feasible states. ^' is evaluated by 
creating a vector of unique states, i/'w for all feasible combinations of J, b, I, F and directly mapping 
them to the values of "f. Define vr as the vector of stationary probabilities where vr {k) is the stationary 
probability of state {k), where A; G {1, • • • , |^|}. Then vr can be evaluated by solving vr = vr^', knowing 
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that vt (i) = 1. 

3) Performance Metrics Calculation: The throughput, R, can be expressed as 

T 



R = W^—^ ^ [PpM-Jk (m) , h (m)) + (1 - pp^a) UbiJk (m) , h (m))], (7) 



k=l m=l 

where Ua (j, i) = I if j £ Fs and i = I, and zero otherwise; Ub (j, = 1 if j € r5 and i = 0, and zero 
otherwise. (m) and (m) are J and / for radio m at state /c, respectively. 
The expected number of colUsions, G, can be expressed as 

G = Y,^''(^) [(1 -Pp.rf) Ua{Jk {m),Ik {m)) + Pp,aUb{Jk {ni),h {m))\. (8) 

k=l m=l 

4) Definitions o/ Pri(-) and Pr2(-).' Pri (•) is defined as in (|4]l. The expression for Pr2(-) is an 
extension to that of the single channel case. Hence, 

M-l 

Pr2 (y) = Prs (^2 (1) 1^1 (M)) Y[ Pts {F2 (m + 1) \F2 (m)), (9) 

m=l 

where Prs(-) is defined as ([5]l replacing /i with x and /2 with y. 

IV. Multi-Stage Spectrum Sensing Algorithm Examples: Analysis 

We present different examples of multi-stages sensing algorithms. As mentioned in Section IIII-All 
Pr3(-) and U{-) are unique to the sensing algorithms and are thus derived here. 

A. Multi-channel Sensing using a Single Narrowband Radio 

Before proceeding with the description of Pr3(-) and U{-), we denote transitions where the SU stays 
on the same channel as Ci = {ci = C2}. Further, denote transitions where the SU switches to the next 

channel as C2 = {ci < M, C2 = ci + 1 or ci = M, 02 = 1}. 

a) No Pre-sensing, Quiet mode (PqQi)' This is the algorithm presented in the IEEE 802.22 stan- 
dard [8|. Since it features a quiet mode but no pre-sensing mode, thus F = {0, • • • , S + 1} and 
Ta = {!,•■■ ; "S* + 1} (note that F5 = {1, . . . , 5} is defined in the same way as in Section IIII-All ). 
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Then 



Prs (3^) 



< 



fl, =0 oris = 0,Ci, 

Pf,s, h{ci) = o,ji e r^, ^2 = ji + i,Ci, 

Pf,t, ^2(ci) = 0, ji = 5+ 1,J2 = 1,(^2, 

1 -Pf,s, hici) = 0, ji € Ts, ^2 = 

1 ^2(ci) = 0, ji =5+ 1,J2 = l,Cl, 

1 -Pm,s, hici) = i,ji e rs,i2 =ii + 

1 -Pm,t, hici) = l,jl = S + 1,J2 = 1,C2, 

^2(ci) = 1, jl G Ts, j2 = l,Cl, 

Pm,t, ^2(ci) = = 5+ 1,J2 = l,Cl, 

Lo, 



(lOa) 
(10b) 
(10c) 
(lOd) 
(lOe) 
(lOf) 
(lOg) 
(lOh) 
(lOi) 
(lOj) 



otherwise. 

We briefly explain the conditions (llOaMlOil ). Transitions to and from the idle mode do not involve 
sensing, hence, Pr3(-) = 1 for transitions where ji = or j2 = as shown in (llOal ). The cases where 
the SU is at any sensing stage and proceeds to the next stage, based on a false alarm or a successful 
detection of a PU, are represented by (llObI ) and (llOfI ). respectively. The cases where the SU is in the 
quiet mode and switches to another channel, based on a false alarm or a successful detection of a PU, 
are represented by (llOcI ) and (lOgi, respectively. The cases where the SU proceeds from an active mode 
of operation to the first stage of sensing, based on correctly detecting the absence of PUs, are represented 
by dlOdl ) and (llOel ). The cases where the SU proceeds from an active stage to the first stage of sensing, 
based on a mis-detection of a PU, are represented by dlOhl ) and (llOil ). respectively. 

Denote transitions where the buffer level remains unchanged and a new frame is produced at the next 
state as C3 = {61 = 625/2 = 1}- Denote transitions where the buffer level decreases by one frame and 



no new frames are produced at the next state as C4 = {62 



1, 61 > 0, /2 = 0}. Then, feasible cases 
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where the SU stays on the same channel, Ci, are 

f JlG{0,--- ,5},i2=0,62 = /2 = 0, 
Jl = S+l,j2 = f2 = B = 0, 
Jl=0,i2 = l,b2 = 0,/2 = 1, 

[/(•) = 1 ^ <! Ji €{!,••• ,5-l},j2 =ii + l,C3 orC4, 

ji = -S, j2 = 5" + 1, /2 = 1, 62 = min(6i + 1, B) 

Jl = S,J2 = 5+ l,/2 = 0,62 = 61,61 >0, 

. 31 e r^, J2 = 1,^3 or C4, 
The feasible case where the SU switches to the next channel, C2, is 



1 



Ji 



S + l,i2 = l,C3orC4. 



(11a) 
(lib) 
(11c) 
(lid) 
(lie) 
(llf) 
(llg) 

(12) 



We explain the conditions behind U{-). A SU that is in the idle mode or is performing multi-stage sensing 
stays at or proceeds to the idle mode, respectively, if no new frames are generated and the buffer is empty, 
as represented in (II lab . A SU with a buffer that is operating in the quiet mode must have at least one 
frame in its buffer, as otherwise, the SU should be in the idle mode. However, if the SU has no buffer, it 
proceeds to the idle mode if there are no new generated frames, as shown in (lllbl) . An idle SU proceeds 
to the first stage of sensing whenever a new frame is generated, the SU transmits the frame and the 
buffer remains empty as represented in (lllcl ). Condition (II Idl) accounts for the case when a SU is in one 
of the first S — 1 sensing stages with a new generated frame and/or has at least one frame in its buffer. 
Then the SU may proceed to the next sensing stage (based on the outcome of the sensing procedure). 
Furthermore, if this SU has a new generated frame, the buffer level remains unchanged; otherwise, the 
buffer level decreases by one. If the SU switches to the quiet mode and has a new generated frame as 



shown in (II let , the buffer level increases by one if the buffer is not full, otherwise, the new frame is 
discarded. If the SU is at the last stage of sensing with at least one frame in its buffer and no new frames 
are generated, the SU may switch to the quiet mode (based on the sensing procedure outcome) and the 
buffer level stays the same as presented in (II Ifb . An active SU with frames to transmit that mis-detects 
a PU or correctly detects the absence of a PU proceeds to the first stage of sensing. In this case, as 
described in (TTgj), if the SU has a new generated frame, the buffer level remains unchanged; otherwise, 
the buffer level decreases by one. If a SU is in the quiet mode and has frames to transmit and the sensing 
procedure results in a false alarm or a successful detection of a PU, the SU switches to a new channel 
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PL' 



and starts at the first stage of sensing. If this SU has a new generated frame, the buffer level remains 

unchanged; otherwise, the buffer level decreases by one. This transition is shown in (fT2l) . 

b) No Pre-sensing, No Quiet Mode (PoQo)' This algorithm aims at reducing the time spent in 
sensing by omitting the pre-sensing and quiet modes of operation. When the outcome of the sensing 
procedure at stage S indicates that a PU is possibly present, the SU switches to a new channel and 
restart the sensing procedure. Then, T = {0, 1, • • • ,5} and Ta = r5 = {1, • • • , S}. Since there are no 
pre-sensing or quiet modes, frames do not get buffered. Accordingly, we assume the absence of a buffer, 
thus, i? = 6i = 62 = 0. Hence 

ji = or j2 = 0,Ci; 

/2(ci) = 0,iiG{l,--- ,5-l},j2=ii + l,Ci or 
hici) = 0,ji = S,j2 = 1,^2; 

^-Pf,s, h{ci) = 0,ii G r5,j2 = l,Ci; 

h{ci) = ,5-l},j2=ji + l,Ci or 

1 Pm,s 5 

hici) = = S,j2 = 1,^2; 
Pm,s, h{ci) = eTs,j2 = l,Ci; 
0, otherwise. 

The derivation of Pr3(-) for FqQo follows that for PqQi. The difference in expressions is caused by the 
absence of the quiet mode. Regarding U{-), feasible transitions where the SU stays on the same channel, 
Ci, are 

fii Gr,i2 = 0,/2 = 0, (14a) 

ii er,j2 = i,/2 = i, (14b) 



(13) 



[/(•) 



1 



< 



bi G {1, • • • , 5 - 1}, j2 = Ji + 1, /2 = 1. (14c) 
The feasible case where the SU switches to the next channel (ci < M, C2 = ci + 1 or ci = M, C2 = 1) is 



[/(■ 



1 



31 = S,j2 = l,/2 = 1- 



(15) 



Since a buffer is unnecessary for this algorithm, a SU with no new generated frames proceeds to the 
idle mode irrespective of its previous state as shown in (I14ab . Condition (I14bl ) indicates that a SU with a 
new generated frame proceeds to the first stage of sensing in case of mis-detecting or correctly detecting 
the absence of a PU. In case of a false alarm or correctly detecting a PU, if the SU is at any of the 
first S — 1 sensing stages, the SU stays on the same channel and proceeds to the next sensing stage as 



17 



presented in (I14cl ). Otherwise, if the SU is at the last stage of sensing, it switches to the next channel 

starting at the first sensing stage as shown in (fTSl) . 

c) Pre-sensing, Quiet Mode (PiQi): As a complement to PoQoj this algorithm features both a pre- 

sensing as well as a quiet mode of operation. If the outcome of the sensing procedure at the pre-sensing 
mode indicates that a PU is possibly present, due to a false alarm with probability pf^f or a sucessfuU 
detection with probability 1 —pm,t, the SU switch to the consecutive channel starting with the pre-sensing 
mode. Otherwise, the SU proceeds to the first stage of sensing at the same channel. In this algorithm, 
r = {0, • • • , S+2} and Ta = {!,••• , S'+2}. The expression for Pr3(-) follows that for FqQi as described 
in ( fTOl ) but, due to the introduction of the pre-sensing mode, the condition on ji in (llOcI ). (llOel ). ( lOgI and 
(II Oil ) changes to ji G {S + l,S + 2}. Also, the condition on j2 in (II Pel ) and ( |10g[ ) changes to j2 = S + 2. 
The expression for U {■) for cases where the SUs stay on the same channel follows that for PqQi as 
given in (fTTb with the following alterations: dllbl ) changes to ji £ {S + I, S + 2}, j2 = f2 = B = 0, and 
(lllcl ) changes to ji = 0,j2 = S + 2, 62 = min(6i + 1,B), f2 = 1. For transitions where the SU switches 
to the next channel 

f ii e {5 + 1, 5 + 2},j2 = S + 2, /2 = 1, 62 = min(6i + 1,B), 
U{-) = 1 ^ < (16) 
[ii e{S + l,5 + 2},i2 = S + 2,/2 = 0, 62 = ^1, bi >o, 

d) Pre-sensing, No Quiet Mode (PiQo): To investigate the effect of introducing the pre-sensing 
mode of operation only, we propose and describe an algorithm that features a pre-sensing mode with no 
quiet mode. Hence, F = {0, • • • ,S,S + 2} and = {1, • • • ,S,S + 2}. The expression for Pr3(-) is as 
given in ( fTOl ) but the condition on ji in dlObl ) and dlOfl ) changes to ji G {1, • • • , S — 1}. The condition 
on ji in (II Pel ) and (II Oil ) changes to ji = S + 2. Finally, the conditions on ji and j2 in (llOcb and (lOg 



change to ji € {S, S + 2} and j2 = S + 2. This is to accommodate for the absence of the quiet mode 
and the presence of the pre-sensing mode. 

Regarding U (•), conditions for transitions for which U {■) = 1 where the SU stays on the same channel 
are given as (ITTI ) with the following alterations. Condition (II Ibb is changed to ji = S + 2, j2 = /2 = 
B = 0. Condition (lllcl ) is altered to ji = 0,j2 = S + 2,62 = min(6i + l,B),f2 = 1. Conditions 
(II lei ) and (II Ifb are dropped as there is no quiet mode. Condition (II Icl ) is altered to accommodate for 
the pre-sensing mode; an idle SU with a new generated frame switches to the pre-sensing mode and 
buffers the new frame if a buffer is available. Feasible transitions where the SU switches to the next 
channel are given as ([T6l ) replacing ji € {S + 1,S + 2} with ji € {S, S + 2}. It implies that if the SU 
is at the last sensing stage or in the pre-sensing mode the SU may switch to the next channel starting 
in the pre-sensing mode depending on the sensing procedure outcome. Since the SU is not allowed to 
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communicate on the channel while in the pre-sensing mode, if a new frame is generated, the frame is 
either buffered, if the buffer is not full, or dropped if the buffer is full or if there is no buffer. Otherwise, 
if no new frames are generated, the buffer level remains unchanged. 

B. Multi-channel Sensing using Parallel Narrowband Radios 

For the parallel narrowband radios case, we describe one algorithm implementation, noting that deriva- 
tions for other algorithms (not presented in this paper) may follow from Section HV-AI Pr3(-), can be 



expressed as 
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pr3w= n p^3,i(3^) n p^3,2(3^) n ^^3,3 (3^), 

mignAi m2&^Q-^^ rrtaefijj 

where Pr3 i(-), Pr3 2(-)> and Pr3 3(-) represent the transition probabilities for radios that are elements of 
sets VLai, f^Qi> and fi/^, introduced in Section HIl-Bli respectively. Pr3^i(-) can be expressed as 



Pr3,i {y) 



1, 

Pf,s: 

1 Pm,si 



P' 



m,si 



Ji {m 

J2 {m 
h {m 

J2 {m 
h {m 

h {m 
h {m 

h {m 



e {!,■■■ ,5-l},J2(m) = 0; 
= 0, Ji (m) G Ts, 

= Ji (m) + 1; 
= 1, Ji (m) G Ts, 

= Ji (m) + 1; 

= 0, Ji (m) G Ts, J2 (m) = 1 or 

= 0, Ji (m) = S, J2 (m) = 0; 
= 1, Ji (m) G Ts, J2 (m) = 1 or 

= 1, Ji (m) = S, J2 {m) = 0; 
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0, 



otherwise. 

The expressions shown in ( fTSl ) follow that for ( fTOl ) where they express whether the active radio proceeds 
to the next sensing stage or back to the first stage of sensing based on the sensing decision and the 
presence or absence of a PU. The only difference is for cases where the radio switches to the idle mode. 
A radio in any of the first 5 — 1 sensing stages can switch to the idle mode in case it is not needed for 
frame transmission. The outcome of the sensing procedure is neglected. For a radio that is in the last 
stage of sensing and is not needed for frame transmission, the radio can switch to the idle mode only 
if the sensing procedure does not detect a PU. Otherwise, the radio switches to the quiet mode. Pr3^2(') 
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can be expressed as 



Pr3,2 (3^) = <^ 



Pf,t, 



I2 (m) — 


u, J2 [m) 


— D i" i, 


I2 {m) = 


1, J2 {III') 


— Q -1- 1 • 


h ijn) = 


0, J2 (m) 


= or 


h (m) = 


0, J2 (m) 


= 1; 


h (m) = 


1, J2 (to) 


= or 


h ijn) = 


1, J2 (to) 


= 1; 


otherwise 
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.0' 

The expressions describe the state transitions for radios in the quiet mode. In case of a false alarm or a 
successful detection of a PU, the radio stays in the quiet mode. Otherwise, the radio proceeds either to 
the idle mode, in case it is not needed for frame transmission, or to the first stage of spectrum sensing. 
Pr3,3(-) can be expressed as 



Pr3,3 C^) 



(20) 



1, J2 (to) = or J2 (to) = 1, 
0, otherwise. 

This expression is for radios that are currently in the idle mode. These radios can either stay in the idle 
mode, if they are not needed for transmission, or proceed to the first stage of spectrum sensing, otherwise. 
No spectrum sensing is done by idle radios. Conditions for feasible state transitions are given by 



1 



Ft, > Ma,M = min{FT, - Ma„B}. 



(21) 



These conditions indicate that the number of active radios should be less than or equal to the number of 
frames available for transmission. Besides, since only Ma, frames can be transmitted. Ft, — Ma, frames 
will need to be buffered. However, as the buffer can only hold B frames, the total number of frames that 
will be buffered equals mm{FT, — Ma,,B}. 



V. Numerical Results 

Due to the vast number of parameters in the proposed model, we focus on the results that capture the 
interplay between PU and SU traffic and sensing algorithms. In specific, we present numerical results 
showing throughput and collision probabiUty as a function of the number of sensing stages. Moreover, 
we investigate the relationship between the sensing time and the performance metrics. Furthermore, we 
demonstrate the relationship between buffer size and the obtained throughput. Finally, we compare the 
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TABLE 11 

Summary of Common Parameters used in Numerical Evaluation 



Descriptive Name 


Parameter Values 


Long-Ts Low-p/ 




= 0.24r, p/,. 


= 0.1 


Short-Ts High-p/ 


Ts 


= o.ir, pf,. 


= 0.36 


Slow PU traffic 


Pp. a — Pptd — 


0.01 


Fast PU traffic 


Pp 


a = 0.5, Pp,d 


= 0.1 



performance of single and parallel narrowband radio architectures. The results for all proposed algorithms 
are shown for single and parallel narrowband radio architectures. 

We adopt a set of parameters and assumptions that are common to all scenarios presented in this section, 
unless otherwise stated. Specifically, we adopt energy detection as the physical layer sensing technique. 
Furthermore, we assume an AWGN channel, where the expressions for the probabilities of false alarm 
and mis-detection for such channel can be found in [19. Eq. (12)] and IHl Eq. (14)], respectivelj]. The 
bandwidth for each PU channel is set to 6 MHz, as in |[T6l Sec. IV], lITOl Sec. IV], the channel throughput 
W = 1Mbps and the SNR at the sensing receiver is set to -10 dB. Slot time is set to T = 1ms, as 
assumed in e.g. |[T8l . Note that we verify all analytical results presented in this section using Monte 
Carlo simulations. 

A. Single Narrowband Radio: Throughput and Collisions versus Number of Stages 

We investigate the effect of varying the number of sensing stages on throughput and collision probability 
for each sensing algorithm in the single narrowband radio architecture case. We assume = 6 channels, 
which represents a typical value for multi-channel networks, cf. G {2,..., 12} |18]. For the SU 
traffic, we assume ps,a = 1 and ps^d = 0, i.e the SU transmits constant bit rate traffic, which saturates the 
channel. We consider SUs without buffer, since our analysis has shown that buffer size does not impact 
throughput when SUs are always transmitting frames. We consider two cases, described below, related 
to sensing to explore the tradeoff between the sensing time and the resulting throughput. We achieve this 
by varying the sensing time and the energy detection sensing threshold while keeping the probability of 
mis-detection constant at pm,s = 0.1. The selected values for pj g are 0.1 and 0.36 which correspond to 

''Note that our model is applicable to any channel model, including the ones discussed in 1101 , 1161 . provided that there is an 
analytical expression linking the sensing time and the probabilities of false alarm and detection for a given SNR. 



21 



Ts = 0.24T and Ts = O.IT, respectively. We denote the former sensing technique by long-T^ low-pj 
and the latter by short-Tg high-pj. The energy detection sensing thresholds for the quiet and pre-sensing 
modes of operation are the same as those for the other stages of sensing. Two realistic PU traffic models 
are considered: (i) one with slow PU traffic (relative to the slot duration) with ^ = pp^d = 0.01, and 
(ii) one with fast PU traffic with pp^a = 0.5 and pp^d = 0.1 lEOl Table 3]. 

Results for the slow PU traffic model are shown in Fig. |l(a)[ Fig. |l(b)[ Fig. 1(c) and Fig. |l(d)[ From 
our framework, it is apparent that the upper bound for R equals ^1 — ^ ^ W . For the slow 

PU traffic model, R is upper bounded by 984.3 kbps. For ideal sensing, that is, if pf^g = Pf,t = Pm,s = 
Pm,t = Ts = 0, R for the four single narrowband radios algorithms approaches the upper bound. Using 
only one stage of ideal sensing results in R that is less than 1% below the upper bound. The loss in 
throughput is caused by the time spent by the SU in the quiet and pre-sensing modes (if they exist 
for the algorithm) and while transmitting on a channel occupied by a PU. However, for single stage 
non-ideal sensing, the yielded throughput for all algorithms is 33-39% below the upper bound for the 
long-Ts \ow-pf sensing technique and 38-53% below the upper bound for the short-T^ high-pj sensing 
technique. 

We now investigate the loss in R for non-ideal sensing for the four algorithms. Considering PqQi and 
PiQo for "S* = 1, it is clear that the quiet mode leads to a higher R compared to the pre-sensing mode. 
This is because for PqQi, a false alarm causes the SU to proceed from sensing stage one to the quiet 
mode. However, the quiet mode features a low probability of false alarm causing the SU to return to 
stage one of sensing, with high probability, and continue communicating on the vacant channel. On the 
other hand, for PiQo, a single false alarm causes the SU to switch to a new channel that might not be 
vacant. PqQo which features no quiet or pre-sensing modes results in a higher R than PiQo for 5 = 1. 
Yet, this comes at the expense of collisions where G for PqQo is 15 and 45 times more than that for 
PiQo for long-Ts low -pf and short-T^ high-pj, respectively. 

Considering multi-stage sensing with 5" > 1 for slow PU traffic, as 5 increases, the impact of false 
alarms on R decreases. This is because more consecutive false alarms have to be generated in order to 
reach the quiet mode (for PqQi and PiQi) or to switch to the next channel (for PqQo and PiQo)- On the 
other hand, increasing S causes the SU to take a longer time, on average, to reach the quiet mode based 
on the successful detection of a PU. This causes G to increase (for most algorithms) and consequently 
leads to a decrease in R. Thus, S* = 2 results in the highest R for most algorithms for long-T^ low-pj. 
However, for short-T^ high-pj, R increases with increasing S due to the dominant effect of false alarms 
on all algorithms. In conclusion, for the slow PU traffic case, multi-stage sensing provides the flexibility 
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Fig. 1. R and G as functions of S for the single narrowband radio case, (a) and (b) sliow R and G, respectively, for slow PU 
traffic and long-Ts low-p / sensing, (c) and (d) show R and G, respectively, for slow PU traffic and short- Ts high-p / sensing, (e) 
and (f) show R and G, respectively, for fast PU traffic and long-Ts low-p/ sensing, and (g) and (h) show R and G, respectively, 
for fast PU traffic and short-Ta high-p/ sensing. The slow and fast PU traffic scenarios and the long-Ts low-p/ and short-Ts 
high-p/ sensing techniques are as explained in Section IV- Al Common parameters: N = 6, T = 1ms, W = 1Mbps, ps,a = 1, 
Ps,d = and B — 0. 
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to decrease the sensing time while increasing the resulting throughput. In our setup decreasing Tg from 
Tg = 0.24r to Ts = O.IT while increasing S to 4 stages leads to increasing R by 14%. Regarding the 
sensing algorithm, the pre-sensing mode results in a considerable decrease in G. Adding the quiet mode 
to the pre-sensing mode leads to an increase in R for small S, and the gain in R is more significant for 
higher probabilities of false alarm. 

Results for the fast PU traffic case are shown in Fig. 1(e) Fig. 1(f) Fig. 1(g) and Fig. |l(h)| As the 



PU traffic becomes faster, the probability that the SU switches to another channel increases. Hence, the 
stationary probability that the SU operates in the pre-sensing or quiet modes increases (for algorithms 
that have these stages) causing the resulting R to decrease. For PqQi, PiQi and PiQo, increasing S 
decreases the stationary probability of operating in the pre-sensing and quiet modes and accordingly, 
causes an increase in R. PqQo exhibits the highest R for the fast PU traffic case, yet, this comes at the 
expense of collisions where G is up to eight times higher than that for algorithms with pre-sensing. For 
PoQo> increasing the number of sensing stages increases the time taken by the SUs to vacate the channel 
for an incoming PU, leading to an increase in G and a decrease in R. Finally, note that the analytical 
results match perfectly with the Monte Carlo simulations. 

B. Single Narrowband Radio: Throughput and Collisions versus Sensing Time Per Stage and Probability 
of False Alarm 

In this section, we analyze the effect of varying Tg on throughput, R, and collisions, G, for the different 
sensing algorithms using the single narrowband radio architecturqj. For this experiment, in contrary to 
Section IV-AI we set 5 = 2 stages, = 3 channels and B = 2 frames to explore the results from the 
perspective of a parameter set that is different from that of the previous section. Furthermore, we set 
Pp,a = 0.01, pp^d = 0.05, ps,a = 0.1 and we set the ratio 1^ (l — ^) ^ ^ (the average generated SU 
traffic) to 500kbps. For the different values of Ts = [50, 500] /is, the energy detection sensing threshold 
is chosen to keep pm,s constant at 0.1, while pf^g = [0.23,0.013]. The unsuccessful frame delivery rate 
is non-zero as frames are lost either because of collisions with PU frames or buffer overflow when there 
are no channel vacancies. For clarity of presentation, we assume a QoS constraint where the maximum 
unsuccessful frame delivery rate is set to 0.1 and we do not present the results that do not satisfy the 
QoS constraint. 



^More results on the effect of Ts on the multi-channel multi-stage spectrum sensing algorithms in the context of energy 
consumption are presented in |1 Sec. IV-B]. 
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Fig. 2. Single narrowband radios case: (a) throughput R and (b) collision probability G as functions of Ts. Common parametrs: 
T = 1ms, 5* = 2, TV = 3, B = 2, Pm,s = 0.1, pp^a ~ 0.01, Pp.d ~ 0.05, Ps,a ~ 0.1 and the generated throughput equals 
500kbps. Note the respective probability of false alarm, Pf,s, for each sensing time are marked at the top of each figure. 



The impact of varying on R is shown in Fig. |2(a)| Results show that for PqQi, PiQi and PiQo, there 
is an optimal value of Ts that maximizes throughput (resulting from the sensing/throughput tradeoff [5], 
0) while for PoQo, R increases monotonically with Tg. This is because increasing Ts decreases the 
probability of false alarm leading to an increase in the utilized spectral vacancy and hence, an increase in 
throughput. On the other hand, increasing Ts increases the number of transmitted frames needed to meet 
the throughput requirement. As a result, for long Ts, increasing Ts increases the probability of buffer 
overflow causing a decrease in throughput. However, this is not the case for PqQo as frames do not get 
buffered, as explained in Section IIV-AI 

The relationship between Ts and G is presented in Fig. |2(b)| Increasing Ts increases the number of 
transmitted frames, to maintain generated throughput at 500 kbps, and thus causes an increase in the 
collision probability. However, shorter Ts results in a higher false alarm probability. This causes the SU 
to switch from a vacant channel to a new channel that might be occupied by a PU resulting in a colUsion. 
The pre-sensing stage decreases collisions caused by switching to an occupied channel. Accordingly, for 
PiQi and PiQo, G increases monotonically with Ts, while on the other hand, the relationship between 
Ts and G is convex and non-monotonic for Ts for PiQi and PiQo- Comparing Fig. |2(a)| and Fig. |2(b)| 
it is clear that the best protocol option (that maximizes throughput, while at the same time minimizes 
collisions) is PiQo- 
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Fig. 3. ii as a function of B for the single narrowband radio case for (a) slow SU traffic and long-Ts low-p/ sensing and (b) 
fast SU traffic and short-Ts high-p/ sensing. The slow and fast SU traffic scenarios are as explained in Section IV-CI and the 
long-Ts low-p/ and short-Ts high-p/ sensing techniques are as explained in Section FV-AI Moreover: S = 2, N — 3, T — 1ms, 
W = I Mbps, pp,a = Pp,d = 0.01. 



C. Single Narrowband Radio: Throughput versus Buffer Size 

For this experiment, we set 5" = 2 stages, = 3 channels, pm,s = 0.1 and pp^a = Pp,d = 0.01. We 
analyze the two scenarios considered eaiiier in Section IV-At (i) one with slow SU traffic (relative to the 
slot duration) with ps^a = Ps.d = 0.01 and the long-T^ low-pf sensing technique, and (ii) one with fast 
PU traffic with ps^a = 0.5 and ps ^ = 0.1 and the short-T^ high-pj sensing technique. 

The impact of increasing S on i? is shown in Fig. [3l Frames are buffered whenever the radio is in the 
quiet or pre-sensing modes. Accordingly, no frames are buffered for PoQo> thus, R is independent of B. 
The buffered frames are transmitted whenever no new SU frames are generated and the SU buffer is not 
empty. For PoQi, PiQo and PiQi, increasing B increases the probability that the new generated frame 
will not be discarded due to the absence of an idle channel. Increasing B increases R where the impact 
of i? on i? is less significant for large values of B. Note, again, the close match between the analysis 
and simulations. 

D. Parallel Narrowband Radios: Throughput and Collisions versus Number of Stages 

We investigate the effect of varying S" on i? and G for parallel narrowband radios. We use the same set 
of parameters as in Section |V-A[ except for = 3 channels. We consider the two PU traffic scenarios 
analyzed earlier: the slow PU traffic with p^.a = Ps,d = 0-01 and the fast PU traffic with „ = 0.5 and 
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Fig. 4. Parallel narrowband radios case: (a) throughput R and (b) collision probbaility G as functions of 5'. Ti and T2 
denote slow and fast PU traffic, respectively, and 5*1 and ^2 denote short-Ts high-py and long-Ts low-p/ sensing techniques, 
respectively. The slow and fast PU traffic scenarios and the long-Ts low-p/ and short-Ts high-p/ sensing techniques are as 
explained in Section FV-AI Moreover: N = 3, T — 1ms, W — 1Mbps, ps,a = 1, Ps.d = 0, B — 0. Note that G exceeds one 
as, in this case, it represents the expected number of collisions for N — 3 per time slot. 



Ps,d = 0.1. We also consider the two sensing options explored earlier: the long-T^ low-pj and the short- 
Tg high-pf sensing techniques. We analyze the scenarios with the four permutations of PU traffic and 
sensing options. The impact of varying 5 on i? and G is shown in Fig. 4(a) and Fig. |4(b)[ respectively. 

Regarding R, increasing S increases R for the four scenarios. This is because more consecutive false 
alarms have to be generated in order for a radio to switch to the quiet mode while the channel is vacant. 
Moreover, the results confirm the intuition that R is directly proportional to N for the tested scenarios. 
Similarly, increasing S increases G for the four scenarios. This is because increasing S increases the 
average time taken by the SU to switch to the quiet mode when a PU arrives. Finally, results clearly show 
the tradeoff between Ts, R and G, and the improvement in R and G caused by employing multi-stage 
sensing. This is demonstrated in the slow PU traffic case with 3 = 4:, where R is higher by 17% and 
G is lower by 35% for scenarios employing the short-T^ high-pj sensing technique compared to those 
employing the long-T^ low-pj sensing technique. As in the previous cases, analysis matches perfectly 
with the simulations. 
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Fig. 5. Performance comparison of single and parallel narrowband radios: (a) successful frame delivery rate and (b) probability 
of collisions per channel per time slot as functions of A'^. Common parameters: T — 1ms, S = 2, ps,a ~ 1, Ps.d ~ 0, B = 0. 
The short-Ts high-p/ sensing technique described in Section IV-AI is adopted; NB: single narrowband radios, WB: parallel 
narrowband radios. 



E. Single versus Parallel Narrowband Radios: Throughput and Collisions versus the Number of Available 
Channels 

Finally, in this section we present the relationship between the performance metrics, R and G, and 
the number of available PU channels, N. At the same time, we compare the performance of single 
and parallel narrowband radios. For this experiment, we set 5 = 2 stages, i? = frames, pm,s = 1 
and pp^a = Pp,d = 0.01. As an example we adopt the short-Ts high-pj sensing technique described in 
Section IV-AI with Tg = O.IT. We consider saturated SU traffic {ps^a = 1 and = 0) for both single 
and parallel narrowband radios. That is, the single narrowband radio generates one frame per time slot, 
as opposed to N frames per time slot for the parallel narrowband radios case. As the throughput of 
the parallel narrowband radios is N times higher than that of the single narrowband radio, the resulting 
number of collisions per time slot is normalized by N to achieve a fair comparison. Moreover, instead 
of using the throughput metric, we introduce a new metric denoted as the successful frame delivery rate, 
which is defined as R/{NW), i.e. the ratio between the successfully transmitted throughput and the 
generated throughput (which in this case equals to the available channel capacity). The number of PU 
channels, N, is increased from two to five and the resulting metrics are presented in Fig. |5] 

The results show that the successful frame delivery rate for all algorithms of single narrowband radios 
is higher than that of parallel narrowband radios. We want to stress, however, that the total (unnormalized) 
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throughput for parallel wideband radios (not shown due to space constraints) is higher than that for single 
narrowband radios. Moreover, the successful frame deUvery rate increases with N for single narrowband 
radios but is constant relative to for parallel narrowband radios. This is because for single narrowband 

radios, increasing the available channels increases the probability of the SU finding a vacant channel 
to occupy. Note that, the successful frame delivery rate is upper bounded by 1 — ( ■ — — I . On the 
Other hand, regarding the parallel narrowband radios case, as we assume saturated traffic, frames are 
generated by the N radios every time slot, hence, the radios can be regarded as operating independently. 
As a result, increasing the number of available channels has no impact on the successful frame delivery 
rate. 

Regarding the normalized probability of collisions per time slot, for parallel narrowband radios it is 
sUghtly less than that for single narrowband radios employing the pre-sensing stage, but significantly less 
than that for single narrowband radios with no pre-sensing. Furthermore, the normalized probability of 
collisions per time slot for parallel narrowband radios do not vary with N due to the independence between 
the parallel radios caused by traffic saturation. Finally, the normalized probability of collisions per time 
slot decreases significantly with N for single narrowband radios with no pre-sensing as increasing the 
available chaimels increases the probability of the SU operating on a vacant chaimel, instead of switching 
to a new chaimel that might be occupied by a PU. 

VI. Conclusions 

We present a unified analytical framework, based on Markov chain analysis, for evaluating multi- 
chaimel multi-stage spectrum sensing algorithms for Opportunistic Spectrum Access networks. Multi- 
stage sensing provides more flexibility in optimizing the systems performance by varying the sensing 
time per stage and the number of sensing stages. Therefore, in the model we consider a variety of 
algorithm design examples that feature prolonged sensing stages before accessing or leaving a channel. 
We also consider the temporal variation in the activity of the primary users as well as the secondary users. 
Furthermore, we analyze operation on multiple chaimels for nodes with single and parallel narrowband 
radios. 

The results demonstrate the tradeoff inherent to multi-stage spectrum sensing between the secondary 
user throughput and the colUsion probabiUty between primary and secondary users. In most cases, 
increasing the number of sensing stages increases both the throughput and the colUsion probability. 
For example, for a single narrowband radio, increasing the number of sensing stages from one to four 
stages, for a slow primary user traffic scenario, results in an increase in throughput of 36%, yet causes an 
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increase in the collision probability by 46%. This tradeoff apphes for both single and parallel narrowband 
radio systems. Moreover, the optimal number of sensing stages has a high dependence on the primary 
user traffic. Generally, having more than two stages does not result in a significant change in throughput 
yet causes an increase in collisions. 

Regarding the sensing algorithms, the pre-sensing mode always results in a decrease in the expected 
number of collisions, while the quiet mode causes an increase in throughput in most scenarios. Results also 
show that increasing the secondary users buffer size results in an increase in throughput, and the increase 
is more significant for fast secondary user traffic. Finally, comparing single and parallel narrowband radio 
systems, single narrowband radios result in a higher successful frame delivery rate as opposed to parallel 
narrowband radios. However, this comes at the expense of the maximum achievable throughput which 
scales with the number of available channels for parallel narrowband radios. 
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